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Diffusion time and splitting of separatrices for nearly 
integrable isochronous Hamiltonian systems 



o 
o 

O I Massimiliano Berti and Philippe BoUe 

Jjz^ • Abstract: We consider the problem of Arnold's diffusion for nearly integrable isochronous Hamilto- 

I nian systems. We prove a shadowing theorem which improves the known estimates for the diffusion time. 

O^ ■ We also justify for three time scales systems that the splitting of the separatrices is correctly predicted 
by the Poincare-Melnikov function. R 
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Riassunto: Tempo di diffusione e splitting delle separatrici per sistemi Hamiltoniani isocroni quasi- 
integrabili. Consideriamo il problema della diffusione di Arnold per sistemi Hamiltoniani isocroni quasi- 
integrabili. Dimostriamo un teorema di shadowing che migliora le stinic sul tempo di diffusione sinora 
note. Giustifichiamo inoltre, per sistemi a tre scale temporali, che lo splitting delle separatrici e corret- 
tamente previsto dalla funzione di Poincare-Melnikov. 



\^ ' 1 Introduction 

We outline in this Note some recent results on Arnold's diffusion for nearly integrable isochronous Hamil- 
CTN , tonian systems: complete proofs and further results are contained in H. We consider Hamiltonians 7i^ 

p^ ' of the form 

P: n^^LU-I+?^ + {cosq-l)+flfiip,q) (1.1) 



with angles variables {}p,q) G T" x T^ and action variables {I,p) £ R" x R^. We assume 
• {HI) There exists 7 > 0, r > n such that \uj ■ k\ > j/\k\'^, Vfc G Z", fc 7^ 0. 



Hamiltonian 7Y^ describes a system of n isochronous harmonic oscillators of frequencies w weakly coupled 
with a pendulum. When fi — Q the energy tuili of each oscillator is a constant of the motion. The problem 
of Arnold's diffusion is whether, for /i 7^ 0, there exists motions whose net effect is to transfer energy 
C^ ' from one oscillator to the others. 

Stemming from Q these kind of results are usually proved noting that, for 11 — Hamiltonian 7i^ 
admits a continuous family of n-dimensional partially hyperbolic invariant tori T/^ — {(ip, I,q,p) € 
T" X R" X T^ X Ri I / = /o, 9 = p = 0} possessing stable and unstable manifolds W(T/J = 14"" (T/ J = 
{{ip,I,q,p) e T" X R" X T^ X Ri I / = /o, p^/2 + {cosq - 1) = 0}. By {HI) aU the unperturbed tori 
7/q, with their stable and unstable manifolds, persist, for fi small enough, being just sligthly deformed. 
The perturbed stable and unstable manifolds W^(7y) and WJ^{Tj^) may split and intersect transversally 
giving rise to a chain of tori connected by heteroclinic orbits. By a shadowing type argument one can 
then prove the existence of an orbit such that the action variables / undergo a variation of 0{1) in a 
certain time T^ called the diffusion time. In order to prove the existence of diffusion orbits following the 
previous mechanism one encounters two different problems 

1) Shadowing theorem; 

2) Splitting of separatrices; 



^Supported by M.U.R.S.T. Variational Methods and Nonlinear Differential Equations. 



By means of a variational technique inspired from H-pl-ipl we give in the next section a general 
shadowing theorem which improves, for isochronous systems, the known estimates on the diffusion time 
obtained in ||^-||ll|-||l^-0 by geometric methods and in [Q by Mather's theory. In section 3, using 
methods introduced in [^ we give some results on the splitting of the separatrices. 

2 The shadowing theorem 

2.1 Perturbation term vanishing along the invariant tori 

We first describe our results when the perturbation term is f{(p, g) = (1 — cosq) f{(p) so that the tori 7/^ 
are still invariant for /i 7^ 0. The equations of motion derived by Hamiltonian Ti.^ are 

(p = uj, /= -/i(I -cosg) 9^/((^), q^p, p ^ sinq - ^smq f{(p). (2.1) 

The angles if evolve as ip{t) — tut + A and then (EJ) can be reduced to the quasi-periodically forced 
pendulum equation 

— q + s\nq = 11 sinq f[ijjt + A), (2.2) 

corresponding to the Lagrangian 

/:^(g,q,i) = ^ + (l-cosg)+M(cos<7-l)/(wi + A). (2.3) 



For each solution q{t) of (2.2) one recovers the dynamics of the actions I{t) by quadratures in (2.1). 

For /x = equation (2.2) possesses the family of heteroclinic solutions qeit) = 4 arctg(exp (i — 9)), 9 G 
R. Using the Contraction Mapping Theorem as in B one can prove that near the unperturbed heteroclinic 
solutions qe{t) there exist, for /i small enough, "pseudo-heterodinic solutions" q'Xg{t) of equation (|2.2|). 



9a 6i(0 ^^^ ^^^^ solutions of (2.2) in each interval (—00, 9) and {9, +00); at time t — 9 such pseudo-solutions 
are glued with continuity at value g^ g{9) = tt and for t —f ±00 are asymptotic to the equilibrium mod 
2tt. Moreover, by a uniqueness property, g^ ^ depends smoothly on {fi,A,9). We can then define the 



function F^ : T" x R ^ R as the action functional of Lagrangian (2.3) evaluated on the "1-bump 
pseudo-heteroclinic solutions" g^g(i), namely 

/9 /'+OC 

^Mieit), <ilAt), t)dt+ C^{q%{t)A%{t). t) dt, (2.4) 

-00 Je 

and the "homoclinic function" G^ : T" — > R as 

G,,{A) ^ F,,{A,0). (2.5) 

There holds F^(A, 9) = G'p(A + Lj9),y9 e R. 

Remark 2.1 The homoclinic function G^ is the difference between the generating functions S^ j (^i 9) 
of the stable and the unstable manifolds ^^'"(7}^) (which in this case are exact Lagrangian manifolds) 
at the section q — tt, namely G^{A) = S~ j {A.tt) — S~^ j (A, tt). Note that G^ is independent of Iq. 

We now give an example of condition on G^ which implies the existence of diffusion orbits. 

Condition 2.1 ("Splitting condition") There exist Aq G T", d>0,0<a<p such that 

• (J) infaBp(Ao) G^c > infBp(Ao) G^ + S; 

• (m) snpB^^Ao) Cm < I + infi3p(Ao) G^; 

• (iii) d{{A e BpiAo) I G^(A) < 6/2+miB,iAo) G^,},{A e Bp(Ao) | G^(A) > 3<5/4+infB^(^„) G^) > 
2a. 



Note that the above "spUtting condition" is clearly satisfied if G^t possesses in Aq £ T" a non- 
degenerate minimum. Moreover Bp{Aq), open ball of radius p, could be replaced by an open subset U of 
T" whose covering set is bounded in R". 

The following shadowing type theorem holds 



Theorem 2.1 Assume (HI) and letG^^ satisfy the splitting condition 2A. Thenyio, Iq withuj-I^ = uj-Iq, 
there is a heteroclinic orbit connecting the invariant tori Tj^ and Tj> . Moreover, there exists C > such 
that Vry > small the "diffusion time" Td needed to go from a rj-neighbourhood ofTig to a r] -neighbourhood 
ofTji is bounded by 

Td < c '^°^^°' pmax (| ln,5|, ^) + C\ Hr^)\. (2.6) 



Remark 2.2 The meaning of (2.t) is the following: the diffusion time Td is estimated by the product of 
the number of heteroclinic transitions k — ( heteroclinic jump / splitting ) — \Iq ~ Io\/S, and of the time 
Tg required for a single transition, that is Td — k ■ Ts- The time for a single transition Tg is bounded by 
the maximum, time between the "ergoditation time" {1/a^), i.e. the time needed for the flow ut to make 
an a-net of the torus, and the time \ lii6\ needed to "shadow" homoclinic orbits for the forced pendulum 
equation. We use here that these homoclinic orbits are exponentially asymptotic to the equilibrium. 

When the frequency vector uj is considered as a constant, independent of any parameter ("a-priori 
unstable" case), it is easy to give a criterion for the splitting condition 2.1 thanks to the first order 
expansion in fx. There holds 

Gf,{A)^ const + fir{A)+0{ji^), (2.7) 

where F : T" ^ R is the Poincare-Melnikov primitive 

r{A) = / (1 " cosqo{t))f{ujt + A) dt. (2.8) 

As a corollary of theorem \2.l\ we have 

Corollary 2.1 Assume (HI) and let T possess a non- degenerate minimum. Then, for ji small enough, 
the statement of theorem ^.l\ holds where the diffusion time is 

Td^O{- log-). (2.9) 

Remark 2.3 The estimate on the diffusion time obtained in /H/ is Td >> 0(cxpl//i) and is improved 
in / pi| / to be Td = 0(exp l//i). Recently in Iw by means of Mather's theory the estimate on the diffusion 
time has been improved to be Td — 0{1/ fi^'^^^). In M/ it is obtained via geometric methods that Td = 
0(l//x'^"'"^). The main reason for which we manage to improve also the estimates o/ fl/ and W is that the 



shadowing orbit of theorem 2.1 can be chosen, at each transition, to approach the homoclinic point Aq, 
only up to the distance a which does not depend of fi. Note moreover that estimate (p^ is independent 
of the number of rotators n. 

Remark 2.4 The above result answers to a question raised in JJ^/ (sec. 7) proving that, at least for 
isochronous systems, it is possible to reach the maximal speed of diffusion /i/|log/i| (moreover ina 
dently on the dimension n). 



2.2 More general perturbation term 

Dealing with more general perturbations /((^s, q) the first step is to prove the persistence of invariant 
tori for /i 7^ small enough. It appears that no more than the standard Implicit Function Theorem is 
required to prove the following well known result. 

Theorem 2.2 Let to satisfy (HI). For fi small enough, for all Iq G R" system TY^ possesses n- 
dimensional invariant tori T^ « 7/„ of the form 

Til = {l = I^ + a^^), 'P = V', <7 = Q^W, P = ^^W, ^ e T"}, 
with Q'^{-), P^{'), a^(') = 0(/i). The dynamics on 7}^ is conjugated to the rotation of speed lo. 

In order to reduce to the previous case we want to put the tori T^ "at the origin" in the {q,p) 
coordinates by a symplectic change of variables. As Tj^ is isotropic, the transformation of coordinates 
(V', J, u, v) -^ (ip, I, q,p) defined, on the covering space R2("+i) of T" x R" x T^ x R^, by 

is symplectic. In the new coordinates each invariant torus 7}^ is simply described by "0 G T", J = 
Iq, u — V = and the new Hamiltonian writes 

{JCf,) JCf, = Ef, + u}- J+ — + {cosu- 1) + Po{fi,u,ij) 

where the perturbation term is 

Po(/i,u» = (cos(Q^ + M)-cosQ'' + (sing'')w + l-cosu)+/i(/(7A,g^ + M)-/(?A,Q^)-9,/(VAg'')w) 

and E^ is the energy of the perturbed invariant torus (a^{^),^p,Q'^{ip), P'^{'il})). Hamiltonian /C^ corre- 
sponds to the quasi-periodically forced pendulum 

— iL + sin u = duPo{fJ',u,ujt + A) (2-10) 

of Lagrangian 

Lf, = — + (l-cosu) -Po{fi,u,ujt + A). (2.11) 



Since the Hamiltonian /C^ is not periodic in the variable u we cannot directly apply theorem 2.1 and the 
arguments of the previous section require some modifications. For ^ small enough there exists, near qg (or 
more exactly its covering orbit in R), a unique pseudo-heteroclinic solution u^ g{t), which satisfies ( 2.10| ) 
in (—00,6') and (0, +00), is glued with continuity at value u^g(6) = n and is asymptotic as t ~* —00 
(resp. +00) to the equilibrium (resp. 27r). Then we define the function T^{A,9) as 



{u- 



J'M-S) = / -^ + {l~cosu%)-Po{ii,u%,ojt + A)dt 



{n%? 



(1 - cosu^ g) - Pi(m, u'X e, wi + A)dt + 27rg^(6'), 



where g^(i) ^ Q^(wi + A) and 

Pi{^,u,ujt + A) = ( cos(q^(t) + m) — cos g^(i) + sin g^(t) (?i — 27r) + 1 — cosu j 

+ m(/M + A, qf^it) +u)- f{ut + A, q^^{t)) - (aj)M + A, qf^it)) (u - 27r) 

We define the "homoclinic function" Gfj.{A) as Gfj.{A) = J-fj,{A, 0). 

The term 27rg^(6') takes into account the fact that the stable and the unstable manifolds of the tori Tj^ 
are not exact Lagrangian manifolds, see ||l5|l . We have 



Theorem 2.3 Assume (HI) and let Gfj. satisfy condition 2.1. Then the statement of theorem \2. \ holds. 
For n small enough we have 

g^,{A)^const + ^lM{A) + 0{^l'^), VAeT" (2.12) 



where M{A) = /+^ f{ujt + A, qo{t)) - f{ujt + A, 0) 



dt. 



Corollary 2.2 Assume (HI) and let M possess a non- degenerate minimum. Then, for fi small enough, 
there exists a diffusion orbit with diffusion time Td = 0(l//ilog(l//i)). 

3 Splitting of separatrices 

3.1 Approximation of the homoclinic function 

If the frequency vector lo — u)^ contains some "fast frequencies" Pi/e'' , & > 0, e being a small parameter, 
the oscillations of the Melnikov function along some directions turn out to be exponentially small with 



respect to e. Hence the development ( p.l2[) will provide a valid measure of the splitting only for /i 
exponentially small with respect to e. In order to justify the dominance of the Poincarc-Melnikov function 
when /i = 0(£P) we need more refined estimates for the error. The classical way to overcome this difficulty 
would be to extend analytically the function Ffj_{A,9) for complex values of the variables, see [|l5[-|10|. 
However it turns out that the function F^ {A, 6) can not be easily analytically extended in a sufficiently 
wide complex strip (roughly speaking, the condition g^ g {Re 9) = t: appearing naturally when we try to 
extend the definition of q^ g to 6* G C breaks analyticity) . We bypass this problem considering the action 
functional evaluated on diS'erent "1-bump pseudo-heteroclinic solutions" Q^ g. This new "reduced action 

functional" Ff^{A,9) has the advantage to have an analytical extension defined for 6 S R + i(— 7r/2,7r/2). 
More precisely let us assume that f{(p, q) = {1 — cosq)f{(p) and that / can be extended to an analytical 
function over D := (R + i[—ai, ai]) x . . . x (R + i[— a„, a„]), for some Oi > 0. Then 



r / 



\k\' 



for aU s e N. 

Define V'o : R ^ R by Vo(i) = cosh^(i)/(l + coshi)^ and set ipe{t) = ipoif - 9). Note that 
/r V'o(i)9o(^) dt ^ 0. By the Contraction Mapping Theorem we find near qe, for /x small enough, 
pseudo-heteroclinic solutions Q^ g{t) and a constant a^ g defined by 



g^ ^j. sinQ'^gf{ujt + A) + a'^gTpg{t) and / (^Q'^g{t) - qe{t)jil:e{t) dt = 0. 






We define the function F^ : T" x R — > R as the action functional of Lagrangian (2^) evaluated on the 
"1-bump pseudo-heteroclinic solutions" Q^g{f), namely 



F^{A,9)= / C,{Q^^^g{t),Q'X,e{t),t) dt 



R 



and G^ : T" ^ R as G^{A) — F^{A, 0). The relation between the functions G^ and G^ is given below 

Theorem 3.1 There exists a smooth diffeomorphism ^'p. '■ T" —* T" of the form ^jf^{A) = A + gf^(A)LU 
with gp{A) : T" — > R satisfying {g^{A), dAfj^iiA)) ^ O as /i — > 0, such that G^ — Gf^ o tp^. 

We now approximate the Fourier coefficients of Gfj,{A) — X^fceZ" Gfc exp*'^'^ with the Fourier coeffi- 
cients of the Poincare-Melnikov primitive T{A) — X^feeZ" ^^^ exp^'^'^. T^ are explicitely given by 

^^•-^^inh(fc.c.f) (^-^^ 



Set||/||=sup^g^|/(A)|. 

Note that the unperturbed separatrix qQ{t) — 4 arctg(exp(i)) can be analytically extended up to 
|/to t\ < 7r/2. Using the Contraction Mapping Theorem, it is possible to extend analytically the function 
Ffj_{A, 0) up to the strip 13 x (R + i(— 7r/2 + 5,tt/2 — (5)), provided /i||/||5~'^ is small. By an estimate of 
F^j^{A, 6) — ^V[A-\-Lo9) over its complex domain and a standard lemma on Fourier coefficients of analytical 
functions we obtain 

Theorem 3.2 There is a constant C such that, for lJ.\\f\\S~^ small enough, Vfc 7^ 0, /c G Z", VS E (0, ^), 

C^2||f||2 



\Gk -f^Tk\< ^jy" exp ( - ^ a,\h\) exp ( - |A; • c.|(| - S)) . (3.2) 

i—l 

3.2 Three time scales 

We consider three time scales Hamiltonians of the form 

L v^ 

n = ^ + e^h + ^ + {coaq - 1) + fi{cosq - l)f{ipi,ip2), h e Tl, h <E IC"-\ n > 2, 
\/£ 2 

namely TL^ with uj^ — (-i=,£°). Such systems have been dealt with for example in |l2| and |17[ |. 

Let /ie~^'^ be small enough. We assume only that / is analytical w.r.t </?2 (more precisely, ai = 
0, and for i > 2, a^ > 0, ai > 7r/2 if a = 0). Set G^{A) = E^^gz Gfei(v42) exp^'^i'^i and T{A) = 
"^k ez -'^'=1(^2) exp**^!"*!. From estimate ( |3.2D we obtain 

Theorem 3.3 For fie^^^^ small there holds 

Gf,{Ai,A2) = Const +(/iro(e,M, ^2) +i?o(e,M, ^2)) +2Re (/iri(e,/x, A2) + i?i(e,M,^2))e'^i 
+ 0(/ie-^/'||/||exp-^) 

where 

Ro{e,fi,A2)=o(fi^\\f\\A and R,{e,fi,A2) = o(^^^^exp-^ 



Remark 3.1 (i) This im prov es the results in \l1{j which require fi = s^ with p > 2 + a. 



(ii) Note that theorem 3^ certainly holds in any dimension, while the results of \12^ , which hold for 
more general syst ems , are proved for 2 rotators only. 
(Hi) Theorem 3.S is not in contradiction with \l3[] . 



This theorem jointly with theorem 3.1 enables us to provide conditions implying the existence of 
diffusion orbits. In fact, if G^ has a proper minimum satisfying condition p.l[ , so has G^. For instance 
we obtain the following result 

Theorem 3.4 Assume that there are A2 and d,c > such that, for all small £ > 0, |ri(A2)| > 
{c/y/e)e~''/^'^^'^ for all \A2 - A2I < d and ro(^2 ± d) > To(A2) + c. Then, for ^e"^/^ small enough, 



condition El is satisfied by G^, with a = Ce^''/^'^^\ 6 == c^/{2y^)e^^/^'^^\ C constant. 



Remark 3.2 In order to prove the splitting of the separatrices using theorem \3.S^ it is necessary, according 
with J7^ and JT^, that 3m, I g Z"~^ such that fo,h fi,7n 7^ 0. 
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